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Hubbard model has been of considerable interests since the discovery of the high temper- 
ature superconductivity In one dimension, the Hubbard model is solvable with Bethe- 
ansatz [Q. The system exhibits an interesting 5*0(4) symmetry [Q. One particular feature 
is that at half-filling, any small positive on-site interaction would make the system a Mott- 
insulator [Q . At less than half filling, the low lying excitations of the system are characterized 
by the universality class of one dimensional interacting electron systems, the Luttinger liq- 
uid. Anderson has suggested that the normal state of the high temperature superconductor 
may be Luttinger- liquid-like, and various studies have been carried out to investigate spin 
and charge dynamics of the Hubbard model and t-J model P,^. 

On the other hand, within the context of problems of interacting electrons, one very 
interesting subject is the theory of itinerant magnetism as a consequence of competition 
of the kinetic energy and the electron-electron interaction, such as the ferromagnetism in 
the one-hole case in the strong interaction limit of the Hubbard model and so on |]6H9|,p!4 



Recently, there are two conjectures on the Hubbard model on the complete graph, based 
on some empirical rules from small size diagonalization. The first was concerned with the 
ground states of the system at less than or equal to half-filling, the second was about the case 
of filling numbers greater than half-filling |TD[. It was conjectured that the ground states of 
the Hubbard model on complete graph are ferromagnetic when the electron number is 
greater than half-filling. In particular, it was conjectured that the ground state is unique 
when Nf, = L + 1. In this work, I provide a proof for this conjecture. Although it is extremely 
simple in my opinions, the details of the proof are still presented carefully, to clarify the 
unclear points existing in previous literatures. 

The Hubbard model defined on a complete graph is given by: 

L 

H = -\t\ cl^Cja + UY^^anih (1) 

o-,l<j^j<-L i=l 

where the size of the lattice is L, the electron number operators are = cJ^Cj^, n^j^ = c||Cj|. 
The electron creation and annihilation operators obey the usual fermion anticommutation 
relations: {cjo-,Cjo-'} = {cia.Cj^i} = 0, {qo-, cj^.,} = 6ij6(j„>. The on-site energy U is assumed 



to be positive, U > 0. For finite U, the free energy of the system can be obtained in the 



thermodynamic hmit L — > oo |11|, however, the wavefunctions of the system were unknown. 



In another hmit where the lattice size is finite and U = oo, the wavefunctions are Gutzwiller 



projected Fermi sea, and they take the Jastrow product form [^,^,0. In the following, 
we are interested in the ground state wavefunctions in the most general case of finite L and 
arbitrary but finite U. 

Assume that the electron number Ne = Yli=iJ2cr=^,icl(jCia = L + 1. Let us work in 
the Hilbert space (denoted by SH ) where the z-component of the total spin operator is 
Sz = 0, (1/2) for even (odd) number of electrons on the lattice. The lowest energy state in 
this space would be the ground state of the Hamiltonian H in the space of iVg = L + 1, 
because the system is rotational invariant. In this subspace SH, one can expand the kinetic 
operator T in terms of its complete eigenvectors: 

T = ^^e„,,,/^|a,s,i^' >< a, (2) 

s,K a 

where the spin quantum number s is summed from to and the quantum number K 

is the total momentum ( T operator is translational invariant ), the Greek letter a stands for 
all other quantum numbers specifying the eigenstates. In this space, the potential operator 
V can also be expanded in terms of a complete basis of vectors spanning the space: 

^ = H ^K-K',o(mod2-K)V{a,s,K;a' ,s,K')\a,s,K >< a ,s,K'\ 

a,cx' s K,K' 

= J2V{A,B)\A><B\>U, (3) 

A,B 

where V{A,B) =< A\V\B >. The operator expansion of V can be transformed from one 
complete basis set to another complete basis set. The last inequality is obtained when using 
the complete basis diagonal in sites. 

Suppose that the state vector (a^, Sg, Kg) is the lowest eigen state of the operator T, with 
the eigenenergy eo = mm{ea,s,K}- If this state vector also satisfies the following property 

V\ag, Sg, Kg >= [/ | , Sg, Kg >, (4) 



then this state is a ground state of the Hamiltonian H = T+V, with eigenenergy Eq = eo+U. 
The operators T and V for the Hubbard model on a complete graph satisfy the condition 
Eq.(|[). It is easy to analyze the spectrum of the T operator in the subspace of fixed total spin 
5* and total momentum K. Introduce the Fourier transformation cl^ = -^Ylm=i^^'^^(4nm 
with m = 0, • • • , M^zi) Then the kinetic operator T becomes 

T = -\t\LY,Eknka + \t\Ne, (5) 

kcr 

where E^ = Skfi and n^o- is the electron number operator in k space. In the Hilbert space 
where Sz = and Sg = {L — l)/2, one can construct the following state vector 

^ 27r(L-l) 

|G>=(E4iC.T)^(4tcSj n 4tIo>, (6) 

1 = 1 fc=27r 
" L 

which satisfies the condition Eq.(^. This wavefunction is thus a ground state of the system 
in the full Hilbert space Ne = L + 1. Obviously, the ground state has eigenenergy 

Eo = -\t\2L + \t\N^ + U, (7) 

where Ne = L + 1. This ground state is ferromagnetic with spin Sg = {Ne — 2)/2, as seen 
from the way it is constructed. The way of construction seems to be the only way that can 
minimize the kinetic energy and interaction respectively. If one tries to minimize H using a 
trial wavefunction which is a linear combination of the lowest energy state of T operator in 
the subspace S = {Ne — 4)/2, the expectation value of H would seem to be greater than Eq, 
because one can have a nonzero chance to have two pairs of electrons on the lattice. This 
can be regarded as a variational argument for the uniqueness of the ground state. 

For the other filling numbers greater than half-filling, similar argument may be carried 
out easily. For instance, suppose that the number of electrons is Ne = L + 2. The ground 
state energy of the system would be 

Eq = -2\t\L + \t\Ne + 2U, (8) 

and the ground state in the Hilbert space of 5"^ = would look like 
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27r(L-l) 

|G',fc'>=(5-)^(44)(cI,^cI,J n 4t|0>, (9) 

where k' can be any value in the range ■ ■ ■ , Hzil^iiLi], xhe ground state wavefunctions 

all have the total spin Sg = {Ne — 4)/2, and ferromagnetism appears in the system. For filling 
numbers Ne = L + f , the ground state energy is Eq = —2\t\L + Ne + fU, and the ground 
state spin is Sg = {L — f)/2. From the above analysis, one can construct the wavefunctions 
explicitly in a similar manner. The rule is pretty straightforward. 

One can generalize the simple argument to discuss ground states of the system defined 
on a complete graph, on which the conduction band interacts with an array of localized 
impurity spins. The hamiltonian is given by 

^ = -1^1 E E cU.-|j|x:?^-5'.(/), (10) 

0-=t,i l<iy^j<L i = l 

where at site i there is a localized impurity spin (spin 1/2). For this impurity lattice 

system defined on the complete graph, the 5*0(4) symmetry no longer holds, which was 
overlooked in previous article |T3|. The electrons interact with the impurity array through 
ferromagnetic exchange. In the Hilbert space where the filling number of the conducting 
electrons is just above half-filling, one can find the ground states of the system. Suppose 
that Ne = J2i=i{Ci']Ci-\ + cJi'^iJ.) = L + 1, one can construct a wavefunction as follows: 

2ir(L-l) 

Ij L — ^ 

|G">=(E4iQT + E^r(/))'^(cWi) n 4tIo>(8) IT, !,■••,!>/, (n) 

i=l i=l k-^ 

which is the ground state of H in the Hilbert space of A^^e = + 1, with eigenenergy given 
by = —\t\2L + Ne\t\ — \ J\{L — l)/4. This ground state is ferromagnetic with total spin 
Sg = iNe + L- 2)/2. 

In summary, in the thermodynamic limit L oo, the free energy of the Hubbard model 
defined on the complete graph consists of two decoupled terms [0], which is an artifact of 
the graph. However, it is still an interesting physical feature that ferromagnetism appears 
in this itinerant electron system, at filling numbers just above the half-filling, because of the 
electron-electron interaction. 
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